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Abstract: The neutron flux at low energy (£■„ < 15 MeV) resulting from the radioactivity of 
the rock in the underground cavern of the India-based Neutrino Observatory is estimated using 
Geant4-based Monte Carlo simulations. The neutron production rate due to the spontaneous fis¬ 
sion of 235,238^^ 232 jj^ («,«) interactions in the rock is determined employing the actual 

rock composition. It is shown that the total flux is equivalent to a finite size cylindrical rock 
{D = L = 140 cm) element. The energy integrated neutron flux thus obtained at the center of the 
underground tunnel is 2.76 (0.47) x 10^^ n cm^^ s^*. The estimated neutron flux is of the same 
order (~ 10^^ n cm^^ s^^) as measured in other underground laboratories. 
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1. Introduction 

Neutrons are known to be an important source of background for experiments like direct dark 
matter searches [1], double beta decay experiments [2], solar neutrino measurements [3], etc. In 
underground laboratories, neutrons originate from the spontaneous fission of U and Th trace el¬ 
ements in the surrounding rock and the (a,n) interactions on the low Z elements present in the 
rock. Neutrons induced by cosmic-ray muon interactions with rock and shielding material are gen¬ 
erally more penetrating because of higher energy {E„ > 15 MeV) but expected flux is ~ 100-1000 
times lower [4,5]. To reach the desired sensitivity in the experiment, the neutron background from 
rock, detector components and cosmic-ray muons should be significantly suppressed. The neutron 
background from rock (E„ <15 MeV) can be reduced by installing passive shield of hydrocarbon 
material or water surrounding the detector. The size and composition of the shield depends on the 
low energy neutron flux. 

In India, a Tin cryogenic bolometer (TIN.TIN) is under development to search for neutrinoless 
double beta decay (Ovj3j3) in ^^^Sn and is proposed to be installed at the upcoming India-based 
Neutrino Observatory (INO) [6,7]. The neutron-induced background {E„ < 18 MeV) in the detector 
and surrounding materials has been studied using neutron activation techniques [8]. To evaluate the 
background levels induced by fast neutrons in the underground cavern and the necessary neutron 
shield requirements, it is essential to have an understanding of the neutron flux in the underground 
cavern. This work describes the estimation of the neutron flux in the cavern from the rock activity 
(spontaneous fission and (a,?i) interacfions). A volume source disfributed uniformly in a finite size 
rock element is considered in the Monte Carlo (MC) simulations, where the strength of the source 
is derived from the INO cavern rock composition. The total neutron flux is then estimated at the 
center of a 12 m long cylindrical tunnel of 2 m radius. Ref. 7 gives the conceptual layout of the 
proposed underground caverns. 
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2. Neutron Yield in the INO Cavern 


The INO cavern will be located in Bodi West Hills (BWH), Madurai, India [9]. The BWH rock 
is mainly Chamockite, which is the hardest known rock, having a density of 2.89 g/cm^. It 
is essential to know the composition of the rock, particularly the content of 235,238jj^ ^^^Th and 
low Z isotopes to determine the neutron production rate. The BWH rock composition is obtained 
from TOF-SIMS (Time of Flight - Secondary Ion Mass Spectrometry) method [10]. Table |T] shows 
various constituent elements of the BWH rock with respective concentrations. Since the SIMS 

Table 1. Elemental distributions of BWH rock obtained with TOF-SIMS method. 


Element 

Concentration 
(% Weight) 

Element 

Concentration 
(% Weight) 

'2C 

2 

39,40,41^ 

4 

23Na 

5 

4«Ca 

9.99 

24,25,26]yjg 

7 

52Cr 

0.49 

22 A1 

25 

56Ee 

1 

28,29,305j 

40 

58,60]^j 

0.55 

31p 

3 

63,65(~;y 

1.2 



'20Sn 

0.6 


method has limited sensitivity >100 ppb, the Inductively Coupled Plasma Mass Spectrometry 
(ICPMS) method [11] was used to reliably extract the 235,238jj^ 232jj^ concentration in the rock. 
It was found that the BWH rock contains 60 ppb of and 224 ppb of ^^^Th. Errors from 
SIMS and ICPMS measurements are about 1%. In the present work, the BWH rock is assumed 
to be a homogeneous mixture of its constituent elements and the 235,238^^ ^^^Th impurities are 
distributed uniformly in the rock. Due to very small isotopic abundance of (0.72%) and low 
SF probability (~ 7 x 10^^ [12]), its contribution to the neutron flux is expected to be negligible. 
Hence, only and ^^^Th decays are considered here. The spectrum of the neutrons emitted in 
Spontaneous fission (SF) is described by an analytic function, known as Watt spectrum [13] and is 
given by: 

W{a,b,E) = Ce^^^‘‘sinh{'/bE) (2.1) 

where C= y/(—)(-^) (2.2) 

The parameters a and b are empirically derived for each isotope [14] (see Table @). 

The (a, n) reaction rate depends on the initial energy of the emitted alpha particle, the reaction 
Q value and the Coulomb barrier. For the present work, the thick target (a, n) reaction yields iN{E) 
- neutrons per MeV) have been taken from Refs. [15,16]. The total thick target neutron yield is 
determined by the sum of the individual element yield weighted by its mass ratio in the BWH 
rock (as per Table |I]). The neutron yields thus obtained, normalized to and ^^^Th content of 
the BWH rock, are shown in Figure [I]. It is evident that the (a,n) component dominates at lower 
neutron energies while the SF dominates at higher energies (E,, >11 MeV). The yield for (a,?i) 
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Table 2. The Watt spectrum parameters for and [14]. 


Isotope 

a 

b 


(MeV) 

(MeV^i) 

238u 

0.64832 

6.81057 

232Th 

0.80 

4.0 



Figure 1. (Color online) Neutron spectra for BWH rock with 60 ppb of and 224 ppb of ^^^Th. EiTors 
are within the point size. 

falls by six orders of magnitude for En ^ 12 MeV and can be neglected. The uncertainties assumed 
in the {a,n) and the spontaneous fission yields are 10% and 8%, respectively [13]. 

3. Monte Carlo Simulations of Neutron Flux 

The underground observatory at the INO site will have a rock cover of at least 1000 m on all sides. 
However, the neutron flux will get attenuated in few meters of rock [17] and hence only a finite 
size rock element will contribute to the neutron flux in the cavern. The size of the rock element 
is selected in such a way that the volume is sufficient to achieve saturation after scattering from 
within the rock element. In other words, neutrons originating outside this volume element have a 
negligible contribution to the flux at surface of the rock element. For simplicity of simulations, a 
cylindrical geometry is considered. As can be seen from Figure |I], neutron source in the rock ma¬ 
terial is very weak and the flux rapidly decreases with energy. This will necessitate the large scale 
simulations and will be prone to errors due to statistical fluctuations. To overcome this problem, 
an alternative approach using MC simulations has been employed. For mono-energetic, uniformly 
distributed, isotropically emitted neutrons within the volume of the rock element, energy spectra 
{N{Ej)) after transmission through a rock element were computed. Since the neutron flux rapidly 


-3- 














decreases with energy (by about ~ six orders of magnitude at £„ = 15 MeV), only <15 MeV 
is considered in the simulations. The transmission factor T{Ej,Ei) - the fraction of neutrons with 
initial energy Ei and emerging with final energy Ej is computed for £■,• = 1 to 15 MeV. The size of 
the rock element was optimized such that T{Ej,Ei) reaches a saturation value for the highest inci¬ 
dent neutron energy £"„ = 15 MeV, over the entire spectrum region (Ej). The neutron flux N{Ej) 
(in units of g^' s^^) originating from the optimized rock element is then computed using: 


N{Ej) = Y,T{Ej,Ei)N{Ei) (3.1) 

i 

where N{Ei) is taken from the total neutron spectrum of the BWH rock given in Figure |T]. The 
neutron flux Ns{Ej) (in units of cm^^ s^^) per unit area at the surface of the detector, over which 
the transmission of neutrons from the rock element is recorded, is given by: 


Ns{Ej) 


PVNjEj) 

A 


(3.2) 


where p is the density of the INO cavern rock, V is the volume of the optimized rock element and A 
is the surface area of the detector kept in front of the rock element. Finally the neutron flux seen by 
the detector is estimated at the center of the tunnel by integrating over the entire cylindrical surface 
of the tunnel. 

The GEANT4-based MC simulation studies have been done using the G4NDL4.2 neutron 
cross-section library [18]. Typically, 10^ events were generated at each incident energy and final 
energy spectra is made with bin widths of 1 MeV. 


4. Neutron Transmission through the INO cavern rock 

To understand the effects of neutron attenuation and scattering as a function of neutron energy 
in the rock, simulations were done for different configurations. For investigating the attenuation 
effect, a mono-energetic point source of neutrons placed behind the center of 40 x 40 cm^ BWH 
rock of different thicknesses was considered. Figure @ shows the neutron transmission probability 
(P(£', )) as a function of rock thickness for £■„ = 1 — 15 MeV. It can be seen that the flux of neutrons 
of £■„ = 15 MeV reduces by about two orders of magnitude after propagation thr'ough 30 cm thick 
BWH rock. It should be pointed out that P{Ei) decreases because neutrons lose energy, and as 
a result the yield of low energy neutrons {N{Ej),Ej < Ei) is enhanced (see Figure However, 
the low energy neutron yield is also affected by scattering in surrounding material. To understand 
the effect of the surrounding material, the rock volume was subdivided into inner cylinder and 
outer cylindrical shell. The inner cylinder was further divided into two cylinders Pi and R 2 with 
diameter D = 30 cm and length L = 30 cm. The outer shell (P 3 ) dimensions were chosen to be 
D = 90 cm and L = 60 cm. These dimensions were chosen on the basis of Figure 0. The schematic 
geometry of these sub-divided volumes (Pi, P 2 ^ P 3 ) is shown in Figure ^ A thin detector of size 
D = 30 cm,L = 0.2 cm is placed on the face of the Pi for recording transmitted neutrons. 
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Figure 2. (Color online) Neutron transmission probability as a function of rock thickness for a mono- 
energetic point source. Lines are only a guide to the eye. 




Figure 3. (Color online) Schematic geometry of rock elements Ri,R 2 and Rt, considered in the simulations. 


Simulations were done for the following three configurations with a uniform volume source : 

• Case 1: 7?i-rock, R 2 and /? 3 -Air : Neutrons in 7?i. This configuration studies the neutron 
propagation through Ri. 

• Case 2: Ri and 7?2-rock, 7?3-Air : Neutrons in 7?2- The case shows the attenuation of neutrons 
after propagation through rock volume . 

• Case 3: R\, R 2 and /? 3 -rock : Neutrons in Ri. The case illustrates the effect of scattering of 
neutrons from adjacent rock volumes R 2 and R^. 

A typical spectrum for the first configuration is shown in Figure ^ 
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Figure 4. Transmission fraction as a function of incident neutron energy Ei — I — 15 MeV and the corre¬ 
sponding transmitted neutron energy {Ej) after propagation through R\ (Case 1). 


The total transmitted neutron spectrum was generated using : 


N{Ej)=Y^T{Ej,Ei)N{Ei) 

i 

where N{Ei) =Nq 5[Ei) (4.1) 

The spectra for all three configurations are shown in the Figure 
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Figure 5. (Color online) (a) Neutron transmission shown for Case 1 (red circle) and Case 2 (blue triangle) 
and (b) Neutron transmission shown for Case 1 (red circle). Case 3 (blue triangle) (see text for details) 


It can be seen from Figure ^a) that the overall contribution from the R 2 is only about 10% 
of that of the R\. Further, reduction in the flux at higher energies (E„ >10 MeV) is greater than 
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that at lower energies. The enhancement in the low energy yield arising due to scattering in sur¬ 
rounding rock material is clearly visible in Figure ^b) in the spectrum corresponding to case 3. It 
is expected that after certain thickness the contribution due to scattering effects will saturate which 
is demonstrated in the next section. 

5. Rock Element Size Optimization 

In order to find the rock thickness where the scattering effects are saturated, simulations were done 
for cylindrical rock elements of different diameters (D - 90, 100, 110, 130, 150 cm). In each case, 
the length of the cylinder was kept same as the diameter to ensure equal transmission length in 
all directions. Since the scattering is expected to enhance the low energy yield, the transmission 
fractions for different energy ranges (0-5, 6-10, 11-15 MeV) were studied for each geometry and 
is shown in Figure 
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Figure 6. (Color online) Transmission fraction of neutron for £„ = 15 MeV as a function of D = L of the 
rock element for the different energy ranges. Each rock element has been scaled to the volume of an element 
with D = L = 90 cm. 


While the transmission fraction is nearly constant for high neutron energies {E„ > 5 MeV), the 
saturation value in the low energies window {E„ = 1—5 MeV) is reached only for rock elements 
size D = L > 130 cm. Hence the optimal size of rock element is chosen to be a cylinder with 
D = L = 140 cm. Figure 0 shows the neutron spectra from this unit rock element used to generate 
transmission matrix T{Ej,Ei) as mentioned earlier. The total neutron flux Ns{Ej) is obtained at 
the surface of this BWH unit rock element using Eq. P?T| . It should be noted that the Ns{E) thus 
obtained represents the flux at any point on the surface of the cylindrical tunnel. 
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Figure 7. Transmission fraction as a function of incident neutron energy Ei — I — 15 MeV and the corre¬ 
sponding transmitted neutron energy {Ej) for rock element D = L= 140 cm. 



Figure 8. (Color online) A schematic diagram of the underground tunnel assumed as a cylinder. The 
surrounding rock is shown by the shaded portion. 


Assuming the experimental hall to be a cylindrical tunnel of radius r = 2 m and L = 12 m (see 
Figure |^, the neutron flux Nf{E) per unit area at its center can be estimated as : 






















































2n 6 


NfiE) =11 N,{E) 

(5.1) 

(p^Oz^-6 


-N,(E)2k I 

(5.2) 

z=-6 


= Ns{E)2% [arctan(^)]^^_6 

(5.3) 

= lMxN,{E) 

(5.4) 


The half length of the cylinder was restricted to 6 m since the contribution to the flux from larger 
lengths will be smaller than 10%. Moreover, additional shield could be provided at this distance in 
the tunnel. Figure ^ shows the estimated neutron flux Nf{E) at the center of the tunnel in the range 
£■„ = 1 — 15 MeV. The energy integrated neutron flux obtained is 2.76 (0.47) x 10^^ n cm^^ s^^ 



Figure 9. (Color online) The estimated neutron flux {Nf{E)) at the center of an underground tunnel at INO 
site. The blue points represent the MC data and the red line corresponds to the fit (see text for details). 


from the rock activity. The errors on the flux include the uncertainties in the MC statistical errors 
and hence the transmission fractions (5%), SIMS and ICPMS measurement and the neutron yields. 
This calculated spectrum can be approximately described by a function: 

N{E)=a exp{-{bE-c)E) (5.5) 

It should be mentioned that neutrons will also originate from the radioactivity of the concrete of 
the experimental hall, the exact composition of which is presently not known. The variation of 
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composition of the rock, if any, is also not considered in the present work. The neutron fluxes from 
the rock measured in other underground laboratories are 3.78 x 10^^ n cm^^ s^^ in the Gran Sasso 
laboratory [13] and 1.72±0.61 (stat.) ±0.38 (syst.) x 10^^ n cm^^ s^^ in the Boulby underground 
laboratory [19] . 

6. Conclusions 

The neutron flux {£„ < 15 MeV) resulting from spontaneous fission and {a,n) interactions for 
BWH rock containing 60 ppb of and 224 ppb of ^^^Th is estimated. Since the neutron source 
in the rock material is very weak and the flux rapidly decreases with energy, an alternative approach 
using MC simulations has been employed. It is shown that only finite rock size contributes to the 
neutron flux at surface and a rock element of size D = L= 140 cm is optimal to evaluate the neutron 
flux. The total neutron flux at the center of a 2 m radius, 12 m long tunnel in the underground cavern 
is determined. The estimated flux at low energy (£■„ < 15 MeV) is 2.76 (0.47) x 10^^ n cm^^ s^^ 
The estimated neutron flux is of the same order (~ 10^® n cm^^ s^^) as measured in other un¬ 
derground laboratories. This will help in the study of neutron-induced background and design of 
neutron shield for TIN.TIN detector. The work will be further extended to study the muon-induced 
neutron flux in the INO cavern. In addition, experiments will be carried out to measure the neutron 
flux when the cavern becomes accessible in future. 
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